Tuning density profiles and mobility of inhomogeneous fluids 
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Density profiles are the most common measure of inhomogeneous structure in confined fluids, but their con- 
nection to transport coefficients is poorly understood. We explore via simulation how tuning particle-wall inter- 
actions to flatten or enhance the particle layering of a model confined fluid impacts its self-diffusivity, viscosity, 
and entropy. Interestingly, interactions that eliminate particle layering significantly reduce confined fluid mo- 
bility, whereas those that enhance layering can have the opposite effect. Excess entropy helps to understand and 
predict these trends. 



Fluids confined to narrow spaces adopt an inhomogeneous 
distribution of density due to the interactions between the fluid 
particles and the boundaries. This density profile provides a 
basic means for characterizing confined fluid structure, and 
in particular how it differs from that of a homogeneous bulk 
system. For simple confined fluids, the density profile can 
be quantitatively predicted using classical density functional 
theory. However, differences between confined and bulk flu- 
ids are not limited to static structure. The former also flow, 
diffuse, and conduct heat at different rates than the latter Un- 
fortunately, a theory that can reliably predict transport coeffi- 
cients has yet to emerge. In fact, even an intuitive understand- 
ing of how the density profile of a confined fluid connects to 
its dynamics is lacking. 

For bulk fluids, semi-empirical structure-property relations 
have helped to correlate and predict transport coefficients (see, 
e.g., idSSl)- Specifically, changes in thermodynamic state 
variables that increase short-range structural order of fluids 
are also known to decrease their mobility in a simple, quan- 
tifiable way. This is true even for systems that exhibit anoma- 
lous dynamical behavior, such as cold liquid water (where vis- 
cosity decreases upon compression) 14|, |5l la] or concentrated 
colloidal suspensions (where interparticle attractions increase 
mobility) 151, 13] ■ Naive extrapolation of this idea might lead 
one to suspect that inhomogeneous fluids with highly struc- 
tured (e.g., layered) density profiles would tend to be more 
viscous and less diffusive than more spatially uniform fluids. 
Is that indeed the case? Here, we explore this issue quantita- 
tively. Specifically, we use molecular simulation to investigate 
the relationship between the transport coefficients of an inho- 
mogeneous fluid and its density profile, the latter of which 
can be modified in a precise way through the interactions of 
the fluid particles with the confining boundaries. 

A key empirical observation motivating this study is the ex- 
istence of an isothermal correlation between the self-diffusion 
coefficient of simple inhomogeneous fluids and excess en- 



tropy (relative to ideal gas), which is approximately obeyed 
across a wide range of confining environments ]|8l Isl llOll . 
Since the magnitude of the excess entropy is itself a mea- 
sure of structural order ll 111 , the aforementioned correla- 
tion is effectively a structure-property relationship. But how 
does excess entropy connect to the density profile? Do fluids 
with more structured density profiles generally have lower or 
higher values of excess entropy when compared to spatially 
uniform fluids? Moreover, can excess entropy be tuned via 
the fluid-boundary interactions to modify the transport coef- 
ficients in a controlled and predictable way? If so, this idea 
might be used to great effect in the engineering of micro- and 
nano-fluidic systems. 

To address these open questions, we turn to the well- 
characterized Weeks-Chandler- Andersen (WCA) model 11211 . 
which is known to capture the entropic packing effects that 
control many properties of dense, atomistic and colloidal 
fluids. The WCA pair potential is defined as (f){r) = 
4e([cr/r]i2 _ [cr/r]^) + eforr < and 0(r) = for r > 

2^/^(7, where r is the interparticle separation. We considerthis 
fluid confined to a thin film geometry between two parallel, 
planar boundaries placed a distance H apart. Particles located 
a distance z from one boundary (0 < z < H) interact with an 
external field 0oxt(z) = 0fw(z) + 4>iw{H ~ z). The single- 
wall potential is given by 0fw(-z) = (2/15)[cr/2;]^ — [(t/z]^ + 
\/l0e/3 H- 0o(^) for z < (2/5)i/V and 0tw(2) = ^0(2) for 
z > (2/5)^/^(7. This represents a WCA 9-3 repulsive bound- 
ary plus an additional term (po{z), which can be used to tune 
the density profile. From here forward, we simplify notation 
by reporting quantities implicitly nondimensionalized by ap- 
propriate combinations of the characteristic length scale a and 
energy scale e (or equivalently k^T, since we set e = k^T for 
all calculations). In the above, /cb is the Boltzmann constant 
and T is temperature. 

Our aim is to investigate, for given values of film thick- 
ness H and average fluid density pavg = p{z)dz. 
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how the details of the density profile p{z) impact relaxation 
processes. This can be done systematically if a suitable set 
of target density profiles can be chosen for study. For fixed 
chemical potential i^l, H, and T, there is a one-to-one mapping 
between 4>o{z) and p{z) (l^. In fact, as explained below, the 
specific (j)o{z) that will produce a given target density profile 
can be determined precisely using Monte Carlo (MC) simula- 
tion methods. Once determined, 4>o{z) can then be imposed 
in an equilibrium molecular dynamics simulation to calculate 
the transport coefficients of the fluid. 

The baseline density profile that we consider is the "nat- 
ural" one for the WCA fluid confined between WCA 9-3 
walls, i.e., the profile adopted by the equilibrium fluid when 
<j>o{z) = 0. It is characterized by a moderately inhomoge- 
neous structure of fluid particles layered parallel to the con- 
fining boundaries. 
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FIG. 1; (a) Natural and flat density profiles p{z), and (b) natural and 
structured density profiles for a confined WCA fluid with average 
density pavg = 0.6 and _ff = 4, as discussed in the text, (c) The 
associated particle-boundary interactions (/)ext(2). 



A second type of profile that we consider is a "flat" den- 
sity distribution, where the layered structure of the confined 
WCA fluid is effectively eliminated by judicious choice of 
(f)Q (z). We define this flat profile to be equivalent to that which 
an equilibrium fluid of non-interacting particles with average 
density pavg would adopt in the presence of the boundary po- 
tential (/)oxt(-z) with 4)o{z) = 0. Fig.fTJa) compares the shapes 
of typical flat and natural density profiles, and Fig.IIlc) shows 
the external potentials 0oxt {z) that produce them for the WCA 
fluid. Interestingly, due to the non-local coupling of ^cxt(^) 
to p{z), the layering structure of WCA particles in the natural 
profile can be generally eliminated by the addition of a non- 
oscillatory, repulsive contribution to the external field. As we 
demonstrate below, flattening the density profile in this way 
has the effect of reducing both the entropy and the excess en- 
tropy of the confined fluid. 

One can imagine that (t)Q{z) could alternatively be chosen 



such that the confined fluid takes on a higher excess entropy 
(and, thus, perhaps a higher mobility) than when i^g = (i.e., 
the natural density profile). To demonstrate this, we investi- 
gate a third class of density distributions which we refer to 
as "structured". The procedure for generating structured den- 
sity profiles is described in detail in the Appendix. In short, 
they correspond to confined fluid states predicted by an ap- 
proximate density functional theory to have particularly high 
values of excess entropy. As is evident in Fig.[TJb), structured 
profiles do have more pronounced layering than their natu- 
ral counterparts, an a posteriori justification for their name. 
Moreover, the MC simulations described below verify that 
they indeed correspond to states with higher excess entropy 
(but lower entropy) than those with the natural profile. 

Here, excess entropy is defined in the usual way as the dif- 
ference between the fluid's entropy and that of an ideal gas 
with the same density profile. To obtain excess entropy data 
for the aforementioned systems, we employed a suite of MC 
techniques. First, we used grand canonical transition matrix 
MC simulations Iil4il with volume V = 1000 to determine 
excess entropy as a function of N (or equivalently Pavg) and 
H for the natural profile case. A detailed explanation of this 
approach can be found in flH- For select values of N and 
H, we then determined the 4)o{z) that generated a target (flat 
or structured) density profile. We subsequently used an ex- 
panded ensemble MC llisll procedure to determine the excess 
entropy of a fluid subjected to this potential. 

To compute i/io (z), we used an efficient nonequilibrium MC 
potential refinement technique recently introduced by Wild- 
ing lfl6ll . The potential was initialized with 0o(^) = and 
subsequently tuned during an eight-stage canonical MC sim- 
ulation. At regular intervals during the ith stage, 4>o{z) was 
incremented by the relative difference between the instan- 
taneous and target density profiles scaled by a modification 
factor Ui = 0.001/2'. Each stage terminated when C, the 
maximum relative difference between the target and aggre- 
gated stage density profiles, dropped below a tolerance of 

C* = 0.01. 

Thermodynamic properties of the confined WCA fluid 
subjected to a nonzero target (t>o{z) were determined through 
canonical expanded ensemble simulations lIlSll . A system 
with a flat or structured profile was related to that with a natu- 
ral profile through a series of subensembles in which the target 
potential 4)o{z) was scaled by a factor A that spanned from 0.0 
to 1.0 in increments of 0.001. A transition matrix MC tech- 
nique similar to that described in UtIi was used to determine 
the relative Helmholtz free energy F(A) of each subensemble. 
The average energy E{\) was also accumulated during a sim- 
ulation. The total entropy difference AS'(A) = S{\) - 5(0) 
between a subensemble characterized by A and a fluid 
with a natural profile was evaluated as A5(A) = 
{[E{\) - E{Q)] - [F{\) - F(0)]). The corresponding 
difference in the excess entropy /S.S°^{\) can be written 
as A5'=-^(A) = A5(A) - A5'^(A). The term A5'^(A) = 
-VH-^ j" [p{z, X) In p{z, X) - p{z,0) In p{z,0)]dz is 
simply the change in entropy of an ideal gas upon changing 
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its density profile from p{z, 0) to p(z, A). 

We computed the transport coefficients of tiie thin films 
via molecular dynamics simulations in the microcanonical en- 
semble using N = 4000 particles and integrating the equa- 
tions of motion with the velocity- Verlet algorithm llSll . A 
time step of 0.0025 was used for simulating the natural and 
flat profile systems, while a shorter time step of 0.0002 was 
employed for the structured profile fluids. Periodic bound- 
ary conditions were employed in the x and y "free" direc- 
tions. We extracted values of self-diffusivity (parallel to the 
walls) D by fitting the long-time (t >> 1) behavior of the 
mean-squared displacement to the Einstein relation for diffu- 
sion (Ar^) ~ 4Dt, where (Ar^) corresponds to the mean- 
squared displacement in the x and y directions. We also calcu- 
lated values of zero-shear viscosity i] using its corresponding 
Einstein relation. 

In order to systematically probe the effects of flattening or 
enhancing the layering of the density profile, we first focus on 
the behavior of the confined fluid for Pavg =^ 0.6 and H ^ 4. 
Let 00, f (2) and 00,3 (-z) represent the contributions to the ex- 
ternal potential that, under these conditions, produce the flat 
and structured profiles, respectively. Starting from the nat- 
ural profile, we incrementally flatten the density distribution 
by choosing 0o(~) = '^f'/'o.f (^) with progressively larger val- 
ues of Af in the range < Af < 1. Similarly, we system- 
atically enhance the layering of the natural profile by setting 
(/)o(z) = As0o,s(-z) with progressively larger values of As in 
the range < As < 1. 

In Fig. Ill we show how these specific ways of modifying 
the density distribution in turn affect the entropy per particle s, 
the excess entropy per particle s°^, the self-diffusivity D, and 
the viscosity rj. As expected, s appears highest for the natural 
profile. This is because the system is virtually athermal when 
00 = due to the steep boundary and interparticle repulsions. 
As a result, its equilibrium (minimum free energy) structure 
also approximately maximizes s relative to other fluid states 
[i.e., other (j>Q{z) and corresponding density profiles] with the 
same pavg and H. The natural profile state does not, however, 
maximize s'^^. Rather, s'^'^ is found to monotonically increase 
with increased "structuring" of the density profile. How can 
we understand these trends? 

The key is to recaU that As''^{X) = As(A) - As°^(A) 
quantifies how changing A and, in turn, the density profile 
modifies the ideal gas entropy, while As°''(A) measures the 
corresponding entropic change associated with the interacting 
fluid's interparticle correlations. With this is mind, it is clear 
that flattening the natural profile will inevitably result in an 
increase in ideal gas entropy [i.e., As"^(Af ) > 0]. Moreover, 
since s generally decreases upon flattening of the natural pro- 
file, it follows that will also generally decrease, reflecting 
an overall strengthening of interparticle correlations. In other 
words, paradoxically, the fluid with the flat density profile ex- 
hibits the highest degree of structural order. 

On the other hand, it is also self-evident that enhancing the 
layering of the natural density profile will decrease the ideal 
gas entropy [i.e., A.s"^(As) < 0]. Since s also decreases in this 
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FIG. 2: Effect of boundary interaction (shape of density profile) on 
(a) excess entropy per particle As™ (A) — s™(A) — s™(0) and en- 
tropy per particle As(A) = s(A) — s(0), (b) self-diffusivity D, and 
(c) viscosity 77 for the confined WCA fluid with pavg = 0.6 and 
H — 4. The centerline corresponds to the fluid with the natural den- 
sity profile of Fig. l(a)and(b). From center to left, the density profile 
is systematically flattened: (poiz) ~ Af0o,f (2), where Af = 1 yields 
the flat profile shown in Fig. 1(a). From center to right, the density 
profile is structured: 0o(z) = Xs(l>o,s{z), where As = 1 produces the 
structured profile shown in Fig. 1(b). Symbols are simulation data, 
and curves are guide to the eye. 



process, the aforementioned entropic penalty can at most be 
partially compensated by an increase in s'^^ due to diminished 
interparticle correlations. In fact. Fig. [2] provides a specific 
example for how "structuring" the density profile can signifi- 
cantly decrease the overall structural order (i.e., increase s'^^) 
of the confined fluid. One should keep in mind that the spe- 
cific structured profiles examined here were selected precisely 
because density functional theory predicted that they would 
have high values of s'^^. Thus, we emphasize that while flat- 
tening the natural density profile can generally be expected to 
decrease s°^, alternative means for embellishing its structure 
at constant Pavg may either increase or decrease s°^. The case 
studied here provides proof of concept for the former. We will 
extensively discuss the latter case in a future publication. 

A second important point of Fig.|2]is that the mobility of the 
confined fluid, as measured by both 77^^ and D, closely tracks 
the behavior of (but not s) for the flattening and struc- 
turing processes described above. A physical explanation is 
that fluid transport parallel to the walls is dominated by inter- 
particle collisions, and hence structural correlations Uj, |2l |3|] . 
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The density profile appears to impact these transport processes 
mostly due to the fact that it modifies the interparticle correla- 
tions, the effect of which is conveniently isolated by comput- 
ing the excess rather than the total entropy of the fluid. 

The above observation that fluids with more uniform den- 
sity profiles can actually have lower and slower dynamics 
than those with strongly layered density profiles (and the same 
Pavg) appears even more general when viewed in the context 
of other recent simulation data. In particular, it is now known 
that the hard-sphere fluid confined between hard walls shows 
lower 0] and slower single-particle dynamics both paral- 
lel 10] and normal lfl9ll to the confining walls when such walls 
are separated by distances that inherently frustrate the ability 
of the fluid to structure into an integer number of layers in its 
density profile. These model predictions of the relationship 
between dynamics and density profiles can now also be read- 
ily tested in experiments, e.g., by using confocal microso py t o 
investigate confined "hard-sphere" colloidal suspensions 12011 . 
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FIG. 3: (a) Excess entropy per particle s'^^, (b) self-diffusivity D, 
and (c) viscosity rj, of the confined WCA fluid versus average density 
Pavg al H — 4. Symbols are simulation data, and curves are guide to 
the eye. 



Finally, to document the generality of the physics discussed 
above, we show in Fig.[3]the behaviors of s°^, rj, and D over 
a broad range of average film densities. Perhaps most striking 
is the observation that, at high density, fluids with structured 
and flat density profiles with the same pavg can differ in both 
D and rj by an order of magnitude. Although the current study 
focuses on equilibrium fluid conditions, the trends evident in 
Fig. [3] suggest that one might even be able to effectively su- 



percool monatomic confined fluids by isothermally modifying 
the external potential in a way that systematically flattens the 
density profile. Finally, one can imagine extending these ideas 
to confined mixtures by using a species-dependent boundary 
potential to tune the partial molar entropies and dynamics of 
individual components. We are currently exploring both of 
these possibilities. 
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APPENDIX 



Protocol for generating high excess entropy density profiles 

Here, we explain how the so-called "structured" target den- 
sity profiles in our study were obtained. As discussed in the 
main text, the idea underlying the structured profile is that 
the corresponding confined equilibrium fluid [with non-zero 
(j)o{z)] should have a high excess entropy relative to the value 
it takes on when it adopts the natural profile [(j)Q{z) = 0] at 
the same H and pavg- Clearly, the name "structured" does not 
follow in any logical sense from the above, but it was rather 
chosen a posteriori to qualitatively reflect the shapes of the 
profiles that we find have this property. 

So how does one search out density profiles associated with 
high excess entropy of the fluid? The approach we adopt here 
is to use classical density functional theory (DFT) for inhomo- 
geneous fluids, which provides a formal, although necessarily 
approximate, connection between excess entropy and the den- 
sity profile. Specifically, we make use of a recent modifica- 
tion 112 in of Rosenfeld's accurate fundamental measure the- 
ory 12211 for confined hard-sphere fluids. Our simulated con- 
fined WCA fluid can be accurately mapped onto a confined 
hard-sphere (HS) fluid with the same effective packing frac- 
tion TrpavgCng/G by assigning the WCA particles an effective, 
HS interparticle diameter <ths- This diameter satisfies the fol- 
lowing "Boltzman factor criterion" 12311 . (f>{r — cths) = !> 
where (f){r) is the WCA interparticle pair potential given in 
the letter 

In the context of DFT, the excess entropy of a confined 
HS fluid can be expressed as, S'^^[p{r)] = —!F'^^[p{r)] = 
— J dr^{r), where <I'(r) is the excess intrinsic Helmholtz 
free energy density. One can further write down an Euler- 
Lagrange equation for maximizing excess entropy, 
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In the above, the second equality follows from the definition 
of the one-body direct correlation function c'^'(r). But how 
can this relationship be understood? 



From the potential distribution theorem H24l 12511 . we have 



c^^(r) = Inpo(r)- where pa{r) is the probability of inserting 
an additional hard-sphere at position r without overlapping 
a particle of the fluid. From this, it is clear that c'-^-' (r) = 
cannot generally be satifisfied for finite Pavg- However, for our 
present purposes, we are not interested in finding states that 
strictly maximize excess entropy. Rather, we want our target 
"structured" profiles to exhibit high values of excess entropy 
relative to that of the natural density profile. Fortunately, for 
that purpose, Eq. lA-ll can still be productively used if recast 
in a different form. 

Specifically, the equation SS^'^[p{r)]/5p{r) = with the 
constraint of constant pavg can be expressed as 



p(z)exp[c(i)(z)] 



(A-2) 



for the slit-pore geometry. Inputting a trial p{z) and the as- 
sociated c(^^(z) predicted from fundamental measure theory 
into the right hand-side of Eq. IA-21 produces a new p{z) on 
the left-hand side. We observed that thus performing Picard 
style iterations on Eq. IA-21 eventually yields density profiles 
that correspond to fluid states with higher excess entropy than 
that of the natural profile. In particular, we carried out these 
iterations with a mixing parameter ll26ll k = 0.005 until the 
relative change in S"^^ (=k~^(S'°'*^„_|_i — S'^^n)) was less than 
0.05. Weighted densities and integrals were calculated using 
the trapezoidal rule on a mesh of O.OOScths- For the highest 
density case, pavg = 0.75, we used the corresponding natural 
density profile as our initial guess. For all lower densities, we 
used the structured density profile at pavg = 0.75 as the initial 
guess. 

All target density profiles and corresponding external po- 
tentials used in this study are available upon request. 
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